EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH (CERN) 




CERN-PH-EP/2012-250 
2012/09/12 



CMS-TOP-1 1-015 



Measurement of the top-quark mass in tt events with 
lepton+jets final states in pp collisions at ^Js = 7TeV 



The CMS Collaboration 



Abstract 



The mass of the top quark is measured using a sample of tt candidate events with 
one electron or muon and at least four jets in the final state, collected by CMS in pp 
collisions at a/s = 7TeV at the LHC. A total of 5174 candidate events is selected from 
data corresponding to an integrated luminosity of 5.0 fb _1 . For each event the mass is 
reconstructed from a kinematic fit of the decay products to a tt hypothesis. The top- 
quark mass is determined simultaneously with the jet energy scale (JES), constrained 
by the known mass of the W boson in qq decays, to be 173.49 ± 0.43 (stat.+JES) ± 
0.98 (syst.) GeV. 
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1 Introduction 

The top-quark mass (m t ) is an essential parameter of the standard model. Its measurement 
also provides an important benchmark for the performance and calibration of the Compact 
Muon Solenoid (CMS) detector HTJ. The top-quark mass has been determined with a high pre- 
cision at the Fermilab Tevatron to be m t = 173.18 ± 0.94 GeV from tt events in pp collisions |2|. 
Measurements have been performed in different top-quark decay channels and using differ- 
ent methods, with the most precise single measurement of m t = 172.85 ± 1.11 GeV being from 
the CDF Collaboration in the lepton+jets channel using a template method J3J . At the CERN 
Large Hadron Collider (LHC), the CMS Collaboration has measured the top-quark mass in the 
dilepton channel 015), where the latest measurement with an integrated luminosity of 5.0 fb _1 
yields m t = 172.5 ± 1.5 GeV 0. The ATLAS Collaboration has published a measurement in the 
lepton+jets channel with an integrated luminosity of 1.04 fb -1 , also using a template method, 
that gives m t = 174.5 ± 2.4 GeV 0. 

In the analysis presented here, we select events containing top-quark pairs where each top 
quark decays weakly via t — > bW, with one W boson decaying into a charged lepton and its 
neutrino, and the other into a quark-antiquark (qq) pair. Hence, the final state consists of a 
lepton, four jets, and an undetected neutrino. The analysis employs a kinematic fit of the decay 
products to a tt hypothesis and two-dimensional (2D) likelihood functions for each event to 
estimate simultaneously both the top-quark mass and the jet energy scale (JES). The invariant 
mass of the two jets associated with the W — > qq decay serves as an additional observable in 
the likelihood functions to estimate the JES directly exploiting the precise knowledge of the 
W-boson mass from previous measurements Q. 

2 The CMS detector 

The central feature of the CMS apparatus is a superconducting solenoid, of 6 m internal diame- 
ter, providing a field of 3.8 T. Within the field volume are a silicon pixel and strip tracker, a crys- 
tal electromagnetic calorimeter (ECAL) and a brass/scintillator hadron calorimeter (HCAL). 
Muons are measured in gas-ionization detectors embedded in the steel return yoke. 

CMS uses a right-handed coordinate system, with the origin at the nominal interaction point, 
the x axis pointing to the center of the LHC ring, the y axis pointing up (perpendicular to the 
plane of the LHC ring), and the z axis along the counterclockwise-beam direction. The polar 
angle, 9, is measured from the positive z axis and the azimuthal angle, </>, is measured in the 
x-y plane. 

The CMS tracker consists of silicon pixel and silicon strip detector modules, covering the pseu- 
dorapidity range < 2.5, where rj = — ln[tan(0/2)]. The ECAL uses lead-tungstate crystals 
as scintillating material and provides coverage for pseudorapidity < 1.5 in the central bar- 
rel region and 1.5 < \rj\ < 3.0 in the forward endcap regions. Preshower detectors are installed 
in front of the endcaps to identify n° mesons. The HCAL consists of a set of sampling calorime- 
ters, that utilize alternating layers of brass as absorber and plastic scintillator as active mate- 
rial. In addition to the barrel and endcap detectors, CMS has extensive forward calorimetry 
that extends the coverage to \rj\ < 5. The muon system includes barrel drift tubes covering the 
pseudorapidity range \rj\ < 1.2, endcap cathode strip chambers (0.9 < \t]\ < 2.5), and resistive 
plate chambers (\r/\ < 1.6). A two-tier trigger system selects the most interesting pp collision 
events for use in physics analyses. A more detailed description of the CMS detector can be 
found in Ref. 111. 
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3 Data sample and event selection 



3 Data sample and event selection 

The full 2011 data sample has been analyzed, corresponding to an integrated luminosity of 
5.0±0.1fb _1 at y/s = 7TeV. Events are required to pass a single-muon trigger or an elec- 
tron+jets trigger. The minimum trigger threshold on the transverse momentum (pj) of an iso- 
lated muon ranges from 17GeV to 24 GeV, depending on the instantaneous pp luminosity. The 
electron+jets trigger requires one isolated electron with pj > 25 GeV and at least three jets with 
p T > 30 GeV. 

We use simulated events to develop and evaluate the analysis method. The tt signal and 
W/ Z+jets background events have been generated with the MadGraph 5.1.1.0 matrix element 
generator [8|, PYTHIA 6.424 parton showering [91 using the Z2 tune IfTUl , and a full simulation 
of the CMS detector based on Geant4 flT]. For the tf signal events, nine different top-quark 
mass values ranging from 161.5 GeV to 184.5 GeV have been assumed. The single-top-quark 
background has been simulated using POWHEG 301 |T2T[T6| , assuming a top-quark mass of 
172.5 GeV. The tt, W/Z+jets, and single-top-quark samples are normalized to the theoretical 
predictions described in Refs. fT7l - [T9"| . The simulation includes effects of additional overlap- 
ping minimum-bias events (pileup) that match their distribution in data. Furthermore, the jet 
energy resolution in simulation is scaled to match the resolution observed in data II20II . 

Events are reconstructed with the particle-flow (PF) algorithm [21J that combines the informa- 
tion from all CMS sub-detectors to identify and reconstruct individual objects produced in the 
pp collision. The reconstructed particles include muons, electrons, photons, charged hadrons, 
and neutral hadrons. Charged particles are required to originate from the primary collision 
vertex, identified as the reconstructed vertex with the largest value of Ep| for its associated 
tracks. The list of charged and neutral PF particles is used as input for jet clustering based on 
the anti-fcj algorithm with a distance parameter of 0.5 p2ll23l| . Particles identified as isolated 
muons and electrons are excluded from the clustering. The momentum of a jet is determined 
from the vector sum of all particle momenta in the jet. From simulation, the reconstructed jet 
momentum is found to be typically within 5-10% of the true jet momentum. Jet energy cor- 
rections are applied to all the jets in data and simulation [20 [. These corrections are defined as 
a function of the transverse momentum density of an event [24-26 1, and also depend on the 
pj and rj of the reconstructed jet. This procedure provides a uniform energy response at the 
particle level with a low dependency on pileup. An additional residual correction, measured 
from the momentum imbalance observed in dijet and photon+jet/Z+jet events, is also applied 
to the jets in data. Finally, the missing transverse momentum is given by the negative vector 
sum of the transverse momenta of all particles found by the PF algorithm. 

Events are selected to have exactly one isolated lepton (£), muon or electron, with pj > 30 GeV 
and 1 7/| < 2.1 and at least four jets with pj > 30 GeV and \rj\ < 2.4. In addition, jets origi- 
nating from bottom (b) quarks are tagged with an algorithm that combines reconstructed sec- 
ondary vertices and track-based lifetime information, the Combined Secondary Vertex Medium 
(CSVM) b tagger described in Ref. H27| . We require at least two b-tagged jets and select 17985 
tt candidate events in data. The estimated selection efficiency for tt signal is 2.3%. From sim- 
ulation, the event composition is expected to be 90% tt, 4% single top quark, 3% W+jets, and 
3% other processes which include the multijet background. Hence, the selection leads to a very 
clean sample of tt events (see also Refs. |[28ll29l '). 
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4 Kinematic fit 

A kinematic fit is employed to check the compatibility of an event with the tt hypothesis and 
thereby improve the resolution of the measured quantities. The fit constrains the event to the 
hypothesis for the production of two heavy particles of equal mass, each one decaying to a W 
boson and a b quark. As indicated above, one of the W bosons decays into a lepton-neutrino 
pair, while the other W boson decays into a quark-antiquark pair. The reconstructed masses of 
the two W bosons are constrained in the fit to 80.4 GeV [7j. A comprehensive description of the 
algorithm and constraints on the fit is available in Ref. |30|. 

The inputs to the fitter are the four-momenta of the lepton and the four leading jets, the missing 
transverse momentum, and their respective resolutions. The two b-tagged jets are candidates 
for the b quarks in the tt hypothesis, while the two untagged jets serve as candidates for the 
light quarks for one of the W-boson decays. This leads to two possible parton-jet assignments 
per event. For simulated tt events, the parton-jet assignments can be classified as correct permu- 
tations (cp), wrong permutations (wp), and unmatched permutations (un), where, in the latter, at 
least one quark from the tt decay is not matched to any of the four selected jets. 

To increase the fraction of correct permutations, we require the goodness-of-fit (gof) probability 
for the kinematic fit with two degrees of freedom P go f = P (% 2 ) = exp (— \x 2 ) to be at least 
0.2. This selects 2906 muon+jets and 2268 electron+jets events for the mass measurement. For 
each event we use all permutations that fulfill this requirement, and weight the permutations 
by their P go f values. In simulation, the fraction of correct permutations improves from 13% to 
44% and the non-tf background is reduced to 4%. 

Figures [T] (a) and [T] (b) show, respectively, the distributions in the reconstructed mass m r ^° of 
the W boson decaying to a qq pair and the mass m\ ec0 of the corresponding top quark for all 
possible permutations before the kinematic fit. The reconstructed W-boson mass m^ c0 and the 
top-quark mass from the kinematic fit mf* , after the P go f selection and weighting, are displayed 
inFigs.[l](c)and[T](d). 

5 Ideogram method 

As the jet energy scale was found to be the leading systematic uncertainty in previous measure- 
ments of m t in this channel, we chose to determine the JES and the top-quark mass simultane- 
ously in a joint likelihood fit to the selected events. The observable used for measuring m t is 
the mass mf* found in the kinematic fit. We take the reconstructed W-boson mass m r ^°, before 
it is constrained by the kinematic fit, as an estimator for measuring in situ a residual JES to be 
applied in addition to the CMS jet energy corrections described in Section |3j This is in contrast 
to a similar measurement from the DO Collaboration [31J, where the kinematic fit is performed 
for different JES hypotheses. 

The ideogram method has been used by the DELPHI Collaboration to measure the W-boson 
mass at the CERN LEP collider |32|. and, at the Fermilab Tevatron collider, by the DO Collabo- 
ration to measure the top-quark mass in the lepton+jets channel |31| and by the CDF Collabo- 
ration to measure the top-quark mass in the all-jet channel [33|. The likelihood in the ideogram 
method is evaluated from analytic expressions obtained and calibrated using simulated events. 
This makes the ideogram method relatively fast and reliable. 

The likelihood used to estimate the top-quark mass and the JES, given the observed data, can 
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5 Ideogram method 



> 

CD 
O 
LO 



CO 

o 

"-4— ' 

CC 

-4—* 

=3 



CD 
Q_ 



CMS, \fs = 7TeV,£+jets 



T 

tt unmatched 
8000- □ « wrong 

tt correct 
7000 - HI tf uncertainty 



| Z+jets 
| | W+jets 
I single top 
• Data (5.0 f b 1 ) 




300 

[GeV] 



> 

CD 

o 

LO 



CO 
£= 

o 

~M — ' 

CO 



CD 
Q_ 



4000 
3500 
3000 
2500 
2000 
1500 
1000 
500 



CMS, \fs = 7TeV,i>+jets 



tt unmatched 
CD tf wrong 
tt correct 
tf uncertainty 



1 I I | I I : 

| Z+jets 
| | W+jets 
I single top 
• Data (5.0 fb 1 ) 




100 



(a) 



200 



(b) 



300 
m reco 



400 

[GeV] 



> 

CD 

<3 3000 

LO 

i2 2500 



CMS, tfs = 7TeV,^+jets 



2> 
'CD 



O 

E 

=5 

0} 



2000 



.2 1500 
cfl 

H — ' 

=3 

E 

CD 
Q. 



1000 



500 



tt unmatched 
I tf wrong 
tf correct 
tt uncertainty 



100 



(c) 



| Z+jets 
] W+jets 
| single top 
Data (5.0 fb 1 ) 



200 



m. 



w 



300 

[GeV] 



> 

CD 

O 

LO 

H — I 

sz 
g> 



CMS, ^s = 7TeV,<?+jets 



CC 

H — I 

Z3 

E 

CD 
Q. 

O 

E 
cn 



1200 r 
1000 7 
800 1 



tt unmatched 
I tf wrong 
tf correct 
tf uncertainty 



| Z+jets 
] W+jets 
| single top 
Data (5.0 fb 1 ) 



.2 600 



400- 



200 




100 



200 



(d) 



300 



400 



ml 11 [GeV] 



Figure 1: Reconstructed masses of (a) the W bosons decaying to qq pairs and (b) the corre- 
sponding top quarks, prior to the kinematic fitting to the tt hypothesis, (c) and (d) show, respec- 
tively, the reconstructed W-boson masses and the fitted top-quark masses after the goodness- 
of-fit selection and the weighting by P go f. The distributions are normalized to the theoretical 
predictions described in Refs. ||T7l - IT9ll . The uncertainty on the predicted tt cross section is indi- 
cated by the hatched area. The top-quark mass assumed in the simulation is 172.5 GeV. 
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be defined as: 

C (m t/ JES|sample) oc C (sample |m t , JES) = fl £ (event|m t JES)™ event 

events 

( \ ^event 
x: c p gof (f)p( m fi>-o| mt jEs) , 
i=l / 

where m t and JES are the parameters to be determined, n denotes the number of permutations 
in each event, and c is a normalization constant. We note that the contribution from background 
is not included in this expression, as the impact of background is found to be negligible after 
implementing the final selections described in Section [4] The ad hoc event weight zv event = 
Ya=\ c Pgoi (0 is introduced to reduce the impact of events without correct permutations. The 
sum of all event weights is normalized by c to the total number of events. 

Due to the mass constraint on the W boson in the fit, mj: 1 ' and are almost uncorrelated and 
the probability P (m^, m^ c °\m t , JES^ for the permutation i factorizes into: 

p(m[;>^°|m t ,JEs) = £/ ; P ; -(mg| mt ,JEs) x P ; (m^|m t JES), 

where fj, with j representing cp, wp or un, is the relative fraction of the three kinds of per- 
mutations. The relative fractions and probability density distributions Pj are determined 
separately for the muon and electron channels from simulated tt events generated for the nine 
different top-quark mass (m t( g en ) values and three different JES values (0.96, 1.00, and 1.04). 
Each of the m^ 1 distributions is fitted either with a Breit-Wigner function, convoluted with a 
Gaussian resolution for correct permutations, or with a Crystal Ball function, for wrong and 
unmatched permutations, for different wz^gen and JES values. The corresponding m^ co distri- 
butions are distorted by the P go f requirement and weighting because permutations with re- 
constructed W-boson masses close to 80.4 GeV are preferred in the kinematic fit. The m^ co 
distributions are therefore fitted with asymmetric Gaussian functions. Figure [2] compares the 
and m^j C0 distributions for the three different kinds of permutations and three choices of 
input top-quark mass and JES to the probability density distributions used for the muon chan- 
nel. A similar behavior for mf 1 ' and is seen for the electron channel. The parameters of all 
fitted functions are parameterized linearly in terms of the generated top-quark mass, JES, and 
the product of the two. 

We obtain separate likelihoods for the muon and electron channels and examine the product of 
these likelihoods to combine the channels. The most likely top-quark mass and JES are obtained 
by minimizing — 21n£ (m tl JESjsample). 



6 Calibration of the ideogram method 

possible biases and for the correct estimation of the statistical uncertainty. For each combina- 
tion of the nine m t/gen values and the three JES scales, we conduct 10 000 pseudo-experiments, 
separately for the muon and electron channels, using simulated tt events. We extract 7ft t , e xtr and 
JES extr from each pseudo-experiment which corresponds to the same integrated luminosity as 
the one analyzed in data. This results in 27 calibration points in the (m t , JES) plane. 

The biases are defined as: 

mass bias = (m t/ext{ — m ti g en ^ , 
JES bias = (JES extr — JES^. 
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6 Calibration of the ideogram method 
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Figure 2: Simulated mf* distributions of (a,d) correct, (b,e) wrong, and (c,f) unmatched tt per- 
mutations, for three generated masses m t/ g en with JES = 1 in (a), (b) and (c), and for three 
jet energy scales with m t/gen = 172.5 GeV in (d), (e) and (f). The vertical dashed line corre- 
sponds to mf 1 = 172.5 GeV. The m]^ co distributions are shown for (g) correct, (h) wrong, and 
(i) unmatched tt permutations for three jet energy scales with m i)gen = 172.5 GeV. The vertical 
dashed lines in (g), (h) and (i) indicate the accepted value of the W-boson mass of 80.4 GeV. All 
distributions are shown for the muon+jets channel. 
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CMS simulation, \s = 7 TeV, n+jets CMS simulation, {s = 7 TeV, e+jets 
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Figure 3: Mean difference between the extracted 7w t/extr and each generated m t/ g en and between 
JES extr and JES gen for the (a) muon channel and (b) electron channel, before the calibration, 
as a function of different generated m tr g en and three values of JES. The colored dashed lines 
correspond to straight line fits, which are used to correct the final likelihoods. The black solid 
line corresponds to an assumption of a constant calibration for all mass and JES points in each 
channel. 



The biases in mass and in JES are shown in Fig. [3] as a function of m t/ g en for the three values of 
JES, and are fitted to a linear dependence to each value of JES. A fit of the calibration points to a 
constant serves as a quality estimator of the overall calibration. Corrections for calibrating the 
top-quark mass m t ca [ and the jet energy scale JES ca j are obtained from the fitted linear functions. 
These final corrections therefore depend linearly on the extracted values of top-quark mass, JES, 
and the product of the two. 

Using pseudo-experiments with the calibrated likelihood, we fit a Gaussian function to the 
distribution of the pulls defined as: 

pull = -. f— , 

where c(m t cal ) is the statistical uncertainty on an individual m t ca [ for a pseudo-experiment 
generated at m tr g en . As depicted in Fig. [4] we find a mass pull width of 1.005 for the muon 
channel and 1.048 for the electron channel Our method slightly underestimates the statistical 
uncertainty of the measurement, and we incorporate the required corrections into the evalua- 
tion of the final likelihoods. 

After applying the single-channel calibration, we again generate pseudo-experiments contain- 
ing both muon and electron events to check the individual calibrations and the calibration for 
the combination of the two channels. As shown in Fig. |5j statistical fluctuations are suppressed 
in the combination, and no additional corrections are needed. 
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6 Calibration of the ideogram method 
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Figure 4: Width of the pull distribution for the calibrated measurement of m t and JES as a 
function of different generated m t/ g en and three values of JES. The muon channel is shown in 
(a) and the electron channel in (b). The black solid lines correspond to fits of constants to all 
calibration points, assuming no dependence on m t or JES. 
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Figure 5: (a) Mean difference between the calibrated and generated values of m t and JES as a 
function of different generated mt,gen and three values of JES for combined lepton+jets events; 
(b) width of the pull distributions for the combined channel after the single-channel calibration. 
The colored dashed lines correspond to straight line fits, the black solid line corresponds to a 
constant fit to all calibration points. The error bars in (a) indicate the statistical uncertainty on 
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7 Systematic uncertainties 

The contributions from the different sources of systematic uncertainties are shown in Table [TJ 
separately for the muon+jets and electron+jets final states, and for the combined fit to the entire 
data set. In general, the absolute value of the largest observed shifts in m t and JES, determined 
from changing the parameters by ±1 standard deviations, are assigned as systematic uncer- 
tainties on the final measurement. The systematic uncertainties considered as relevant for this 
measurement, and the methods used to evaluate them are described below. 



Fit calibration: We propagate the statistical uncertainty on the calibration to the final measured 
quantities. 

b-JES: The difference in the jet energy responses for jets originating from light (uds) or bottom 
quarks or gluons studied in simulation indicate that the response to b-jets lies between 
the jet responses to light-quark and gluons |20|. Hence, the uncertainty assumed for the 
flavor dependence of the JES determination, which covers the transition from a gluon- 
dominated to a light quark-dominated sample, also covers the difference between light 
quarks and bottom quarks. Thus, all the momenta of all b jets in simulation are scaled up 
and down by this flavor uncertainty. 

pj- and f/-dependent JES: As we measure a constant jet energy scale we have to take into ac- 
count the influence of the pj- and //-dependent jet energy uncertainties. This is done 
by scaling the energies of all jets up and down according to their individual uncertain- 
ties [20] . We take the largest difference in the measured top-quark mass and JES (com- 
pared to the expected average JES shift of 1.6%) as a systematic uncertainty. 

Lepton energy scale: We shift the muon and electron energies in simulation up and down ac- 
cording to their respective uncertainties. 

Missing transverse momentum: In addition to propagating the jet and lepton energy scale 
uncertainties, we scale the unclustered energy up and down by 10%. 

Jet energy resolution: The jet energy resolution in simulation is degraded by 7 to 20% depend- 
ing on tj to match the resolutions measured in data in Ref. [20]. To account for the reso- 
lution uncertainty, the jet energy resolution in the simulation is modified by ±1 standard 
deviations with respect to the degraded resolution. 

b tagging: The nominal requirement on the CSVM tagger is varied in order to reflect the un- 
certainty of the b-tag efficiency of 3% |I27| . 

Pileup: The effect of the pileup events is evaluated by weighting the simulation to provide as 
many additional minimum bias events as expected from the inelastic pp cross section. To 
cover the uncertainties associated with the number of pileup events, the average number 
of expected pileup events (9.3 for the analyzed data) is changed by ±5%. 

Non-tt background: After final selection, background fractions of 1% W+jets and 3% single 
top-quark events are expected from simulation. We conduct pseudo-experiments with 
2% W + bb and 6% single top-quark events and take the difference to the result without 
background as a systematic uncertainty. 

Parton distribution functions: The simulated events are generated using the CTEQ 6.6L par- 
ton distribution functions (PDFs) |34). The uncertainty on this set of PDFs is given in 
terms of variations of 22 orthogonal parameters that can be used to reweight the events, 
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7 Systematic uncertainties 



Table 1: List of systematic uncertainties for the muon+jets and electron+jets final states, and for 
the combined fit to the entire data set 
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0.004 


0.23 


0.004 


b tagging 


0.11 


0.001 


0.15 


0.001 


0.12 


0.001 


Pileup 


0.07 


0.002 


0.08 


0.001 


0.07 


0.001 


Non-tt background 


0.10 


0.001 


0.16 


0.000 


0.13 


0.001 


Parton distribution functions 


0.07 


0.001 


0.07 


0.001 


0.07 


0.001 


Renormalization and 


0.23 


0.004 


0.41 


0.005 


0.24 


0.004 


factorization scales 


ME-PS matching threshold 


0.17 


0.000 


0.15 


0.001 


0.18 


0.001 


Underlying event 


0.26 


0.002 


0.24 


0.001 


0.15 


0.002 


Color reconnection effects 


0.66 


0.004 


0.39 


0.003 


0.54 


0.004 


Total 


1.06 


0.008 


1.00 


0.007 


0.98 


0.008 



resulting in 22 pairs of additional PDFs. Half of the difference in top-quark mass and JES 
of each pair is quoted as a systematic uncertainty. 

Renormalization and factorization scales: The dependence of the result on the renormaliza- 
tion and factorization scales used in the tt simulation is studied by changing the nomi- 
nal renormalization and factorization scales simultaneously by factors of 0.5 and 2. This 
change in the parameters also reflects the uncertainty on the amount of initial- and final- 
state radiation. 

ME-PS matching threshold: In the tt simulation, the matching thresholds used for interfacing 
the matrix-elements (ME) generated with MadGraph and the pythia parton showering 
(PS) are changed from the default value of 20 GeV down to 10 GeV and up to 40 GeV. 

Underlying event: Non-pertubative QCD effects are taken into account by tuning PYTHIA to 
measurements of the underlying event [10]. The uncertainties are estimated by compar- 
ing two tunes with increased and decreased underlying event activity relative to a central 
tune (the Perugia 2011 tune compared to the Perugia 2011 mpiHi and Perugia 2011 Teva- 
tron tunes described in Ref. II35II ). 

Color reconnection effects: The uncertainties that arise from ambiguities in modeling color 
reconnection effects [36 1 are estimated by comparing in simulation an underlying event 
tune including color reconnection to a tune without it (the Perugia 2011 and Perugia 2011 
noCR tunes described in Ref. II35I ). 

Differences in the systematic uncertainties calculated for the muon+jets and electron+jets final 
states are consistent with arising from statistical fluctuations due to the limited sizes of the 
studied samples. The total systematic uncertainty is dominated by effects that cannot be com- 
pensated through a simultaneous determination of m t and JES. Besides the JES for b quarks, 
the uncertainty in color reconnection dominates. For the sample without color reconnection 
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Figure 6: (a) The 2D likelihood (— 2Alog (£)) measured for the £+jets final state. The ellipses 
correspond to statistical uncertainties on m t and JES of one, two, and three standard devia- 
tions, (b) The statistical uncertainty distribution obtained from 10 000 pseudo-experiments is 
compared to the uncertainty of the measurement in data of 0.43 GeV. 



effects, the mean of the m^ co distribution shifts relative to the reference sample, leading to an 
additional uncertainty on for the simultaneous fit. 

8 Measurement of the mass of the top quark 

Using the selected samples, we measure: 

^+jets: w t =173.22 ± 0.56 (stat.+JES) ± 1.06 (syst.) GeV, JES = 0.999 ± 0.005 (stat.) ± 0.008 (syst.), 
e+jets: m t =173.72 ± 0.66 (stat.+JES) ± 1.00 (syst.) GeV, JES = 0.989 ± 0.005 (stat.) ± 0.007 (syst.). 

The combined fit to the 5174 £+jets events in the two channels yields: 

mt = 173.49 ± 0.43 (stat.+JES) ± 0.98 (syst.) GeV, 
JES = 0.994 ± 0.003 (stat.) ± 0.008 (syst.). 

The overall uncertainty of the presented measurement is 1.07 GeV on the top-quark mass from 
adding the components in quadrature. The measured JES confirms that obtained from events 
with Z bosons and photons |20|. 

Figure [6] (a) shows the 2D likelihood obtained from data. As depicted in Fig.|6](b), the uncer- 
tainty of the measurement agrees with the expected precision from the pseudo-experiments. 
As the top-quark mass and the JES are measured simultaneously, the statistical uncertainty on 
tn t combines the statistical uncertainty arising from both components of the measurement. 

We estimate the impact of the simultaneous fit of the jet energy scale by fixing the JES to unity. 
This yields m t = 172.97 ± 0.27 (stat.) ± 1.44 (syst.) GeV. The larger systematic uncertainty stems 
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8 Measurement of the mass of the top quark 



CMS, is = 7TeV,f+jets 



> 




Figure 7: Distribution of the reconstructed top-quark mass after the kinematic fit, for the per- 
mutation with the lowest x 2 1 in the cross-check analysis. The distributions are normalized 
to the number of events observed in data. The top-quark mass assumed in the simulation 
is 172.5 GeV. The rightmost bin also contains the overflow. The inset shows the cross-check 
likelihood measured for the £+jets final state. 

from a JES uncertainty of 1.33 GeV and demonstrates the gain from the simultaneous fit to m t 
and JES. 

As a cross-check of the event selection and the mass extraction technique, a second analy- 
sis is performed. The mass extraction technique used in this analysis is very similar to the 
CMS measurement of the mass difference between top and antitop quarks l|29l . The events are 
required to pass a lepton+jets trigger and fulfill the same lepton and jet requirements imple- 
mented in the main analysis, except for a lower threshold on the muon transverse momentum 
of px > 20 GeV. In addition, the requirement on the number of b-tagged jets is lowered to at 
least one. The same kinematic fit is employed as for the main analysis. All possible permuta- 
tions of the four jets of largest pj that have a fit x 2 < 20, corresponding to a goodness-of-fit 
probability of P go f = 4.5 x 10~ 5 , are accepted, yielding 54899 selected events. This fit to m t 
employs just the standard jet energy corrections, equivalent to setting JES=1. New likelihood 
functions are formed that take account of the contributions from background resulting from 
less stringent selection criteria. In addition, the permutations are weighted with the probabili- 
ties for b tagging p9|. Figure [7] shows the mass mf 1 after the kinematic fit for the permutation 
with the smallest x 2 i n each event and the likelihood obtained from data. After applying the 
calibration, we measure a top-quark mass of m t = 172.72 ± 0.18 (stat.) ± 1.49 (syst.) GeV, which 
is consistent with the result of the main analysis but has a larger uncertainty as expected. 

We use the BLUE method 113 711 to combine the result presented in this letter with the measure- 
ments in the dilepton channel in 2010 [4j and 2011 |5|. Most of the systematic uncertainties 
listed in Table [l] are assumed to be fully correlated among the three input measurements. Ex- 
ceptions are the uncertainties on pileup, for which we assign full correlation between the 2011 
analyses but no correlation with the 2010 analysis, since the pileup conditions and their treat- 
ments differ. In addition, the mass calibration, the statistical uncertainty on the in situ fit for 
the JES, and the data-based background normalization in the dilepton analyses are treated 
as uncorrelated systematic uncertainties. As a cross-check, the correlation coefficients for the 
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correlated sources are varied simultaneously between 1 and 0, changes of < 0.14 GeV are ob- 
served in the combined result. The combination of the three measurements yields a mass of 
m t = 173.32 ± 0.27 (stat.) ± 1.02 (syst.) GeV. It has a % 2 of 1.05 for two degrees of freedom, 
which corresponds to a probability of 59%. 



9 Summary 

The complete kinematic properties of each event are reconstructed through a constrained fit 
to a tt hypothesis. For each selected event, a likelihood is calculated as a function of as- 
sumed values of the top-quark mass and the jet energy scale, taking into account all possible 
jet assignments. From a data sample corresponding to an integrated luminosity of 5.0 fb -1 , 
5174 candidate events are selected and the top-quark mass is measured to be m t = 173.49 ± 
0.43 (stat.+JES) ± 0.98 (syst.) GeV. This result is consistent with the Tevatron average 0, and 
constitutes the most precise single measurement to date of the mass of the top quark. 
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